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ABSTRACT 

Aims. We study the possible detection of and properties of very high-energy (VHE) y-ray emission (in the energy band above 
100 GeV) from high redshift sources. 

Methods. We report on the detection of VHE y-ray flux from blazars with redshifts z > 0.5. We use the data of Fermi telescope in 
the energy band above 100 GeV and identify significant sources via cross-correlation of arrival directions of individual VHE y-ray s 
with the positions of known Fermi sources. 

Results. There are thirteen high-redshift sources detected in the VHE band by Fermi/LAT telescope. The present statistics of the 
Fermi signal from these sources is too low for a sensible study of the eff'ects of suppression of the VHE flux by pair production 
through interactions with Extragalactic Background Light photons. We find that the detection of these sources with ground-based 
y-ray telescopes would be challenging. However, several sources including BE Lacs PKS 0426-380 at z = 1.11, KUV 00311-1938 
at z = 0.61, B3 1307+433 at z = 0.69, PG 1246+586 at z = 0.84, Ton 116 at z = 1.065 as well as a flat-spectrum radio quasar 4C 
+55.17 at z = 0.89 should be detectable by HESS-II, MAGIC-II and CTA. A high- statistics study of a much larger number of VHE 
y-ray sources at cosmological distances would be possible with the proposed high-altitude Cherenkov telescope 5@5. 
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1. Introduction 

y-ray emission from distant blazars is suppressed by the ef- 
fect of pair production through interactions of these y-rays with 
the low-energy photons forming the Extragalactic Background 
Light (EBL) (|Gould & Schreder, 1967; Kneisk reTaLl [2004 
'Stecker et al., 2006; Mazin & Raue, 2007; Franceschini et al. 



,20Q8^^Gilmore et al.^ 2012^ ). This prevents observations of high- 
redshift sources using the technique of imaging the Cherenkov 
emission produced through y-ray induced air showers, used 
by the ground-based Cherenkov telescopes HESS, MAGIC and 
VERITAS ( Aha ronian et a l. , 2008 ). Most of the VHE y-ray loud 
blazars, detected by these telescopes, are situated in the local 
Universe, at (several) hundred megaparsec distances, in the red- 
shift range z ~ - 0.2[] Only one source at redshift z > 0.5, 
3C 279, was possi bly detected by the MAGIC telescope during 
a flaring activity ( [Albert et al.[ [20Q8| ). Another source at red- 
shift 0.444, 3C 66A, was detected by VERITAS ( [Acciari et aL| 
[2009] ). One more relatively high redshift source, PG 1553+113 
(at z > 0.4, reported by'Danfor th et al.[p 01Q )) has been detected 
by MAGIC and VERITAS ( Albert et al., 2006 ; 0^[2QTT ). 

Measurement of the effect of suppression of the VHE y- 
ray flux from low-redshift blazars is commonly used for the 
estimation of the density of the EBL in the local Universe 
([Kneiskeetal.|[2 QQ4; Stecker et al. , 2006; M azin & Raue| |2007| 
[Franceschini et al.^, 2008^ Gilmore et al.^ 2012| ). Constraints on 
the EBL density and spectral characteristics, extracted from the 
VHE blazar observations, are useful for understanding the cos- 
mological evolution of galaxies. Observations of the suppressed 



VHE y-ray emission from blazars at non-negligible redshifts can 
provide a measurement or constraint not only of the EBL density 
in the local Universe, but also new information on the cosmolog- 
ical evolution of the EBL, which is largely uncertain. 

Observations of high-redshift sources in the VHE band 
would also provide valuable information on the cosmological 
evolution of the VHE y-ray loud blazars (and, more generally, 
y-ray emitting active galactic nuclei, AGN), which is also highly 
uncertain. This information is important, in particular, for an 
und erstanding of the origin of Extragalacti c y-ray background 
( Ab do et al.[ |2Q10at [Neronov & Semikoz| [T012). Most of the 
VHE y-ray emitting blazars at zero redshift belong to a High- 
energy peaked BL Lac (HBL) sub-class of the blazar population. 
Negative cosmological evolution of this sub-class was reported 
based on the observations in the visible and X-ray band (IGionmii 



et al. 1999 ). Such puzzling evolution is, apparently, opposite to 



the general positive evolution of blazars and radio galaxies (BL 
Lac parent AGN population) with t he cumulative power of the 
sources incre asing as (1 + ^ > ( (Hodge et al. 2009 Sadler et 



lal., 2007; Smolcic et al.l 



2009 ). Independent verification of this 



http : //tevcat . uchicago . edu 



hypothesis using y-ray observations would be possible only if a 
significant amount of VHE emitting blazars could be detected in 
the redshift range z ~ 0.5 - 1. 

Taking into account the importance of the study of the VHE 
y-ray emission from distant AGN, we report here on the detec- 
tion of VHE y-ray signals from sources at redshifts z > 0.5 
by the Large Area Tel escope (EAT) on board of Fermi satellite 
( Atwood et al. 2009 1). The eff'ective area of LAT, about 1 m^, 
is several orders of magnitude smaller than that of the ground- 
based y-ray telescopes, so that LAT has detected only one- 
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to- several photons from the brightest extragalactic VHE y-ray 
sources in ~ 4 years of observations. However, an extremely low 
level of the background (including the residual cosmic ray back- 
ground and Galactic / extragalactic y-ray background) makes 
even a several-photon signal in the energy band above 100 GeV 
significant. This fact has been used by [Neronov et al.| ( |2010[ 



2011 ) for the search of extragalactic VHE y-ray sources via the 
analysis of clustering of VHE y-rays on the sky, based on the 
method first proposed by jOorbunov et al.| ( |2005 |). In this paper 
we use the correlation of arrival directions of VHE y-rays de- 
tected by Fermi with the positions of high-redshift blazars to 
identify high-redshift sources of VHE y-rays and to study their 
spectral characteristics. 



2. Data selection and data analysis 

For our analysis we have used Fermi pass 7 data within the time 
window from August 4, 2008 up to June 11, 2012. We have con- 
sidered two sub-classes of the "ULTRACLEAN" events ("ev- 
class=4" in the photon selection routines): the event from the 
classes 65311 and 32543. Those classes have the best angular 
resolution in Fermi, corresponding to a "superclean class" in 
pass 6 data for £^ > 100 GeV. Using the method of Gorbunov 
et^(|2005) we correlated the arrival directions of photons with 
£" > 100 GeV with positions of the sources listed in the two- 



year Fermi catalog ( [Nolan et"ST| |2012| ) with high Galactic lati- 
tude \b\> 10°. This analysis is si milar to one done for the IFGL 
catalog by [Neronov et al.| (20ir), where a catalog of 75 objects 
with £" > 100 GeV was obtained. 

There are 1036 front and 888 back converted photons at high 
Galactic latitude \b\ > 10° and E > 100 GeV in the selected data. 

We have measured the PSF for the selected events via the 
correlation of photon arrival directions with the sky positions of 
known gamma-ray blazars. The size of the 68% containment cir- 
cle for the event classes 6531 1 and 32543 is 0.08° for front pho- 
tons and 0.21° for back photons with energies above 100 GeV. 
In order to keep the analysis methods consistent with those used 
by|Neronov et al.| ( |2011| ), we use 0.1° and 0.2° distance bins for 
correlation analysis using front and back photons. 

We took all high Galactic latitude Fermi LAT sources with 
known redshifts and selected high redshift sources with z > 0.5. 
The original two-year Fermi catalog (Nolan et al.| |2012^ does 
not contain information on the redshift of the sources. Instead, 
we have co nsidered the redshifts listed i n the 13th Veron catalog 
of BL Lac s (|Veron-Cetty & Veron 2010) and in the Rome blazar 
catalog of Massaro et al. (|2007). We cross-checked the informa- 
tion on the source redshifts with the Fermi AGN catalog ( Abdo 
[et a l., 2010b) and with the information given in the NECj^and 
SIMBAE[^ databases. In this way out of 414 Fermi BL Lacs we 
selected 55 with z > 0.5. 

For the 55 selected BL Lacs we found six photons correlat- 
ing with six objects within 0.1 degree, while a total background 
of 0.06 was expected. The probability that this signal occurred 
by chance was P ~ 6 • 10"^^ The six sources contributing to the 
correlation are listed in the upper part of the Table [T] Within 0.2 
degrees from the source positions, (but outside the 0.1 degree re- 
gion of the sources) the expected background was 0.14 photons, 
while 4 photons were observed from 3 additional BL Lacs. 

One of these three additional sources was KUV 0031 1-1938, 
with two associated photons. The chance coincidence probabil- 
ity for the two photons to arrive within 0.2 degrees from the 



source position was 8 x 10"^, which corresponds to a 4.5(T de- 
tection of the source above 100 GeV. The chance coincidence 
probability for the remaining two photons to arrive within 0.2° 
from two other BL Lacs was P ~ 10"^. We list these sources 
in the lower part of the Table [T] (two last lines), but one should 
remember that they might, in principle, be "fake detections" due 
to the chance coincidence of the background photon arrival di- 
rections. 

From the 351 Flat Spectrum Radio Quasars (FSRQ) we se- 
lected 269 objects with z > 0.5. Within 0.1 degree from these 
objects we found four photons from four different sources, while 
a total background of 0.4 was expected. The probability that this 
happened by chance was P ~ 8 • 10""^. Again, these objects were 
given lower priority in Table 1, similar to that for the case of BL 
Lacs with photons only associated within 0.2 degrees. We did 
not consider correlations with FRSQ objects within larger an- 
gles because of the high level of chance coincident background. 

The complete list of sources with redshift z > 0.5 which 
correlate with the arrival directions of the VHE y-rays in LAT is 
shown in Table [T] 

For the spectral and timing analysis presented in Figs. 1-13, 
we have used the Fermi Science Tools We have calculated the 
source spectra using two complementary techniques: the likeli- 
hood analysis and the aperture photometry methods. The two 
methods give consistent results. The likelihood method being 
more reliable at low energies (0.1 - 10 GeV) where the LAT 
instrument has a relatively wide point spread function. The aper- 
ture photometry method enables to properly take into account 
the Poissonian statistics of the signal at low photon count rates 
in the 10-300 GeV range. 



3. Properties of individual high-redshift sources 

In this section we study in more detail the spectral and timing 
properties of the high-redshift VHE y-ray sources listed in Table 
[T] to understand if some of the lower significance sources in 
Table [T] among the FSRQs and the two BL Lacs listed in the last 
two rows of the table could possibly be random coincidences of 
arrival directions of the VHE y-rays with the source positions on 
the sky. We also investigate the possibility for the detection of 
these sources using ground-based y-ray telescopes. 



3.1. Firm detections 
3.1.1. TXS 0138-097 

The redshift of this so urce has been obta ined through spectro- 
scopic measurements ( Stocke et al. 1997| . At the same time, the 
SIMBAD database entry for the source gives a redshift of 1.03, 
based on the data fro m the SDSS photometric catalog ( tSDSS| 
Data release 7 |2009|). The detec tion of absorption lines at red- 
shift z = 0.5 (Stick el et al.[p^9T| for this source imposes a lower 
limit on its redshift, so that in any case the source belongs to the 
z > 0.5 source sample in which we are interested in this paper. 

The upper panel of Fig. [T] shows the lightcurve of the source 
in the energy band above 1 GeV. From this lightcurve, one can 
see that the source behaviour is stable on the time scale of al- 
most four years of Fermi exposure. The lower panel of the same 
figure shows the source spectrum in the 0.3-300 GeV energy 
range. The spectrum is well described by a simple power law of 
the form dN/dE ~ E~^ with spectral slope F ^ 2, suppressed at 



^ http://ned.ipac.caltech.edu/ 

^ http://simbad.u-strasbg . fr/simbad/ 



http ://f ermi .gsfc.nasa. go v/s sc/data/analy sis/ 
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Table 1. List of high-redshift blazars emitting in the energy band £^ > 100 GeV. Columns A^o.i and give the numbers of VHE 
photons within an angular distance of 0.1° and 0.2° from the source position. P is the chance coincidence probability for the VHE 
photons to be found within the 0.1° or 0.2° circles from the source labels "f ' and "b" mark photons pair converted in the front- 
and back- layers of the LAT detector. Column A^so-ioo gives the number of photons associated to the source in the lower energy 
band 30 - 100 GeV. Emax is the maximal energy of photon associated to the source. L/Lurk 421 is the source luminosity in units of 
Mrk 421 luminosity at ^ > 100 GeV 



high-energies through absorption on the EEL. The model spec- 
trum shown in Fig.[T]is calculated assuming a source redshift of 



3.1.2. PKS 0426-380 



0.733 and using the EEL model of ,Franceschini et al.|pQQ8 



The 100-200 GeV band source flux is consistent with the expec- 
tation based on the Franceschini et al. model. 



Of course, the statistics of the Fermi signal is largely insuffi- 
cient for a sensible study of the details of the shape of the spec- 
trum in the 100 GeV band, which is determined by the effects 
of propagation through the EEL. Such a study would be possi- 
ble only with a ground-based y-ray telescope, which provides a 
much larger collection area for y-rays and, consequently, much 
higher source signal statistics. The lower panel of Fig. [T] shows 
a comparison of the measured source spectrum with the sensi- 
tivity limits of different existing and future ground-based y-ray 
telescopes. Unfortunately, the model source flux is below the 
sensitivity limit of the HESS-II telescope, which will start op- 
eration in 2012. The estimate of the HESS-II sensitivity shown 
in this plot has been taken from |Masbou| pOlO^) . The sensitiv- 
ity curve found in Masbou| (|2010|) is given in terms of the inte- 
gral flux. Conversion of the sensitivity into the differential flux 
sensitivity depends on the slope of the source spectrum. Fig. [T] 
shows three difl'erent differential sensitivity curves correspond- 
ing to three different values of the slope of the source spectrum. 
One can see that within the energy range covered by HESS-II, 
the spectrum is steep, with a slope close to F = 4, so that the 
source flux is below the HESS-II sensitivity level such that the 
source is not expected to be detectable. 



The expected performance of the next-generation ground- 
based y-ray telescope CTA is not better than that of HESS-II 
for energies below 100 GeV. Because of this, the prospects for 
the source detection with CTA are also not promising. Contrary 
to this, however, the sensitivity of a ground-based y-ray tele- 
scope optim ized for the 10 GeV energy band, 5@5 ( Aha ronian| 
et al. 2Q00| ) is sufficient for the source detection in the energy 
range above 10 GeV and up to the sharp EEL-induced cut-off' at 
100 GeV. The diff'erence in the performance of 5@5 and CTA in 
the 10-100 GeV energy band turns out to be crucially important 
for the possibility of the study of VHE y-ray emission from this 
high-redshift source. 



This bright EL Lac is situated at a redshift of z = 1.11 (|Heidt 
et al. 2004 ). The source is strongly variable in the GeV energy 
band, see upper panel of Fig. [2] Several active episodes, during 
which the flux increased by up to an order of magnitude, have 
been observed during the four years of Fermi's operation. 

The source spectrum averaged over the whole observation 
period extends as a soft powerlaw with F > 2 well into the 
100 GeV energy band at the level of ~ 10"^^ erg/cm^s, see lower 
panel of Fig. |2] At this flux level, the source should be readily 
detectable in the 30-100 GeV energy range with the HESS-II 
telescope and, in the future, with CTA. 

We have verified that the source flux did increase by a fac- 
tor of ~ 2-3 during the flaring activities, without any sig- 
nificant change in its spectral shape. This implies that the de- 
tectability of the source from the ground should be facilitated if 
the next ground-based observations are triggered during the next 
increased activity episode (which could be traced based on the 
Fermi monitoring of the source). 

3.1.3. B2 0912+29 

The SIMEAD database ascribes a redshift z = 1.521 to the 



source, with a reference to the SDSS photometric survey (SDSS 
"Data release 7 2009| ), but stating that the redshift is spectro- 
scopic. At the same time, the NED database does not ascribe a 
redshift to the source. 

The source spectrum can be described by a relatively hard 
powerlaw with slope F ^ 2 up to ~ 200 GeV energy, showing no 
sign of suppression by the EEL. Such an unabsorbed powerlaw 
shape of the spectrum would pose a problem for modeling of the 
source if it is indeed at a redshift of z = 1.5. At this redshift, 
the suppression of the flux due to absorption on the EEL should 
start already below 100 GeV, so that the flux measurements in 
both 50-100 GeV and 100-200 GeV bands would be in tension 
with the model, see Fig.|3] 

The bottom panel of Fig. |3] shows a comparison of t he sensi- 
tivi ty of th e MAGIC telescope in stereo operation mode ( Colin et 
[aLl 120091 with the level of the source flux for different assumed 
values of the slope of the powerlaw spectrum F. One can see that 
if the source flux is indeed at the level measured by Fermi, rather 
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Fig. 1. Lightcurve (top) and spectrum (bottom) of TXS 0138- 
097 at redshift z = 0.733. Black curve in the bottom panel 
sho ws a powerlaw type model spectrum absorbed on the EBL 
fromlFranceschi ni et al.| p008). Also shown in the bottom panel 
are sensitivity curves o f different ground based y-ray te lescopes: 
5@5 ( [Aharonian et a l., 2000), HES S-I an d HESS-II ( |Masboul 
|2QTq1 ) and CTA (|CTA ConsortiunJ [20TTJ. Black and red data 
points in the lower panel show the spectrum extracted using 
the likelihood analysis and aperture photometry method, respec- 
tively. 



than the theoretically expected flux from a source at the redshift 
1.521 suppressed by the absorption on EBL, this source should 
be detectable by MAGIC. 
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Fig. 2. Lightcurve (top) and spectrum (bottom) of PKS0426-380 
at the redshift z= 1.11. Notations are the same as in Fig.[T] 



for such a situation. A proper study of the shape of the spectrum 
and of the efl'ect of the EBL on it could be done only with a ded- 
icated instrument with a very low energy threshold, like 5@5. 
From Fig. |4] one can see that the 5@5 instrument would be able 
to detect the source starting from an energy of ^ 10 GeV. 



3.1.5. PG 1246-H586 



3.1.4. Ton 116 

The redshift of the source cited in the NED database is z = 1 .065, 



with a reference to the SDSS measurement ( [SPSS Data release 
2007|). A t the same time, the Fermi AGN catalog ( |Abdo "et 
2010b| ) list this source as having an unknown redshift. The 
source exhibits a hard spectrum with slope F ^ 1.7 (see Fig. 
|4]), so that in spite of the high redshift of the source, the flux 
in the 100 GeV energy range should be at the level detectable 
by MAGIC and, in the future, by CTA. The source spectrum is 
expected to become very steep just above 100 GeV, so that its 
signal in MAGIC and CTA would be dominated by the contri- 
bution from the lowest energy bin just at the threshold energy 
of the instrument. It is not clear what quality of measurement of 
the shape of the spectrum in the VHE band could be achieved 



The source redshift listed in both the NED an d SIMBAD 
databases is z = 0.847 with a reference to SDSS ( |SDSS Data] 
release 6[|2007| ). In the Fermi AGN catalog the source redshift is 



considered as unknown ( j Abdo et^L} |20 1 Ob| ) . The source exhibits 
a steady behavior, remaining stable over the four-years of Fermi 
monitoring (see the top panel of Fig. [5]) and the source spectrum 
is hard with a slope slightly harder than F = 2 up to energies 
above 100 GeV, where the spectrum starts to be afl'ected by the 
efl'ect of suppression on the EBL (bottom panel of Fig. [5|. The 
sensitivity of the MAGIC telescope would be sufficient for the 
detection of the source, as would the sensitivity of CTA. Similar 
to the case of Ton 116, a proper study of the shape of the VHE 
y-ray spectrum could be done only with a dedicated low energy 
threshold instrument, like 5@5, which would be able to detect 
the source starting from an energy -10 GeV. 
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E [eV] 

Fig. 3. Lightcurve (top) and spectrum of B209 12+29 at the sup- 
posed redshift z = 1.521. Notations are the same as in Fig. [T] 
Blue curves show the sensitivity of the MAGIC telescope in the 
stereo observation mode, from [Colin et al.| ( |2QQ9| ). 



E [eV] 

Fig. 4. Lightcurve (top) and spectrum of TON 1 16 at the redshift 
z = 1.065. Notations are the same as in Figs. [T] and [3] 



3.1.6. B3 1307+433 

The redshift quoted for this source in t he NED database i s 
z = 0.69, based on the SDSS measurement ( [Plotkin et aLl|2008| ). 
At the same time, SIMBAD q uotes a redshift of z = 2.159, also 
based on the SDSS data ( SDSS Data release 7||25Q9|). The spec- 
trum of the source, shown in Fig.|6] is a hard (F < 2) powerlaw 
up to ~ 200 GeV. If the source were at a redshift of 2.159, as 
mentioned by SIMBAD, the detection above 100 GeV would be 
in contradiction with the expected level of suppression of the 
source flux due to the interactions of y-rays with the EBL pho- 
tons. At a redshift of 0.69, however, the Fermi measurement is 
consistent with the theoretical expectation (see Fig. [6]). 

Verification of the source spectral properties in the 100 GeV 
band, therefore, would allow clarification of the issue of the 
source redshift via observations with the existing ground based 
telescopes (MAGIC and Veritas). However, high quality mea- 
surements of the shape of the spectrum would require a next 
generation instrument, like CTA and 5@5. 



3.1.7. KUV 00311-1938 

As is mentioned in the previous section, this source is individ- 
ually detected above 100 GeV with significance 4.5cr. The red- 
shift of this BL Lac is z = 0.61 ( [Giommi et al.[|20()5] ). Its spec- 
trum is characterized by a hard powerlaw type with F < 2 from 



0.3 up to ~ 200 GeV. The VHE data point in the 100-200 GeV 
band matches exactly that expected from an extrapolation of the 
low-energy spectrum, with no signature of the EBL suppression 
noticable at this energy at the Fermi data quality level. The Fermi 
measurement, however, is consistent with the assumption of a 
powerlaw suppressed by the absorption on EBL. The EBL eflfect 
is expected to appear mostly above 200 GeV, in the energy band 
readily accessible by the HESS telescope. The source flux above 
100 GeV is sufficiently high for the detection with HESS-II and, 
possibly, HESS (see Fig.[7|. 



3.2. Possible detections 

In this subsection we provide the details on sources for which 
the correlation of the arrival directions of £^ > 100 GeV pho- 
tons with their source positions is less significant. We notice that 
some "chance coincidence" sources can possibly be singled out 
based on their spectral characteristics (the extrapolation of the 
low-energy spectrum to the VHE band is not consistent with the 
source flux level in the VHE band implied by the detection of 
the VHE photon). However, this criterium alone can not be used 
to firmly reject the source due to the possible existence of pecu- 
liar source spectra (e.g. due to the presence of diff'erent emission 
components in diff'erent energy bands). 
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Fig. 5. Lightcurve (top) and spectrum of PG 1246+586 at the Fig. 6. Lightcurve (top) and spectrum (bottom) of KUV 00311 



redshift z = 0.847. Notations are the same as in Figs. [T] and |3] 



3.2.1. 4C +55.17 

4C +55.17 is a source from the FSRQ subclass of blazars, i.e. 
from the sub-class which gives only a minor contribution to the 
set of known VHE y-ray sources in the local Universe. The ap- 
pearance of this source, along with the three other FSRQs in 
the z > 0.5 VHE source sample is perhaps important since it 
may provide an indication of the higher FSRQ contribution to 
the VHE y-ray source population in the cosmological past. 
The source's redshift is z 



0.8955 ( [Adelman-McCarthy et| 
aL||2008| ). The source spectrum is softer than F = 2 for GeV en- 
ergies and above, as expected from FSRQs which have on aver- 
age softer spectra than BL Lacs (see Fig. [8]). An extrapolation of 
the £^ < 100 GeV band spectrum into the VHE band is consistent 
with the estimate of the VHE flux of the source based on the de- 
tected VHE photon. This indicates that the VHE y-ray emission 
from the source is real. However, detection of the source from 
the ground would be a challenging task for the existing ground- 
based telescopes like MAGIC (see Fig. [8|. At the same time, 
the source should be readily detectable by the next-generation 
instruments CTA and 5@5 (Fig. [8]). 

3.2.2. TXS 1720+102 

This is another representative of the FSRQ sub-population in the 
high-z VHE source sample. The source is at a redshift z = 0.732 
( [Afanasiev et al.[ |2005| ). Its spectrum is softer than F = 2 in 
the energy band below 100 GeV. The overall source flux level is 



at the redshift z = 0.69. Notations are the same as in Figs. [T] and 



somewhat low compared to e.g. 4C +55.17, so that the source 
would be hardly detectable by MAGIC, Veritas and/or CTA (see 
Fig. |9]), unless the spectrum has a "bump"-like feature in the 
100 GeV band. The quality of Fermi data is not sufficient to 
judge whether this is the case or the Fermi detection in the E > 
100 GeV is due to the up-fluctuation of the photon count rate. 
The £" > 100 GeV data are well consistent with the extrapolation 
of the lower energy powerlaw attenuated the EBL (see Fig. [9]). 

3.2.3. PKS 1958-179 



The spectral characteristics of this FSRQ at a redshift ofz = 0.65 
( [Barkhouse & Hall 2WT) appear somewhat peculiar with a soft 
spectrum below an energy ~ 30 GeV and a possibly a new hard 
component appearing above the 30 GeV band, see Fig.fTOl The 
statistics of the Fermi data is not sufficient to draw definitive 
conclusions about the existence of two separate components in 
the spectrum. In principle, a single powerlaw with a photon in- 
dex of about F = 2 provides a satisfactory fit to the data above 
an energy of 1 GeV. Within such a spectral model, the detection 
of the source in the energy band above 100 GeV is consistent 
with the assumption that the powerlaw extends up to the VHE 
band, if the absorption of the VHE y-rays through interactions 
with the EBL photons is taken into account. The source could be 
marginally detectable with the ground based y-ray telescopes, 
including both existing and future facilities (see Fig.fTOl). 
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redshift z = 0.61. Notations are the same as in Fig.[T] 



the redshift z = 0.8955. Notations are the same as in Figs.fTland 



3.2.4. PKS 2142-75 



This FSRQ at z = 1 . 139 ( [Jauncey et al. 1978 ) also has a peculiar 
time-average spectrum in the 0.3-300 GeV energy band, with a 
soft component dominating at energies below 30 GeV and pos- 
sibly a new hard component appearing above 30 GeV, see Fig. 
1 1 If not for this hard component, the source would not be de- 
tectable at energies above 100 GeV. During the period of Fermi 
observations, the source has exhibited several pronounced flares 
(Fig. 1 1 upper panel). The complicated shape of the spectrum 
might be due to the fact that difl'erent components contribute to 
the source flux during the flare and quiescent periods. However, 
with only two photons above 30 GeV, one could not make any 
firm conclusions on the reality and possible nature of the "hard 
excess". It is well possible that the VHE photon coming from 
within 0.1° of the source is a background photon falling coin- 
cident to the source (as discussed in the previous section, one 
out of the four FSRQs listed in Table [T] could possibly be a false 
detection in the VHE band). 

3.2.5. RGB J0250+172 



Contrary to the sof t spectra FSRQs, RGB J0250+172 at z = 1.1 
( [Bauer et aLl|2000| ) is characterized by a hard y-ray spectrum and 
is a good candidate for being a real VHE y-ray source. However, 
due to the lower overall flux normalization, the source should 
only be marginally detectable by both Fermi and by the current 
and next-generation ground-based y-ray telescopes (see Fig. [12]). 



Unexpectedly, a photon with energy 358 GeV is detected from 
the source direction. This detection is somewhat puzzling. On 
the one side, the estimate of the source flux in the 200-400 GeV 
band, based on this photon, is consistent with the powerlaw ex- 
trapolation of the lower-energy spectrum, without an account of 
the EBL-induced suppression of the source flux (see Fig. [12]). On 
the other side, if the EBL suppression is taken into account, the 
measured source flux is orders of magnitude above the expected 
flux from a source at redshift 1.1. Resolution of this inconsis- 
tency might be (a) that the VHE photon from the source direction 
is a background photons, with the chance probability 0.8%, (b) 
that the energy of this particular photon is over-estimated by the 
LAT data analysis software or (c) that the redshift of the source 
has not been correctly measured. 

This contradiction could be easily resolved via observation 
of the source with existing ground-based y-ray telescope(s) like 
MAGIC. If the source flux in the 200-400 GeV band is indeed at 
the level of ~ 10"^^ erg/(cm^s), it should be readily measurable 
by MAGIC. If the flux is at the level expected based on the esti- 
mate of the model calculation taking into account the EBL sup- 
pression, the source is likely not detectable neither by MAGIC 
nor by CTA. It would be, however, detectable by 5@5 due to the 
lower energy threshold of the instrument. 
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Fig. 9. Lightcurve (top) and spectrum (bottom) of TXS 
1720+102 at the redshift z = 0.732. Notations are the same as in 
Figs.[T]and[3l 



3.2.6. PKS 1130+008 

The redshift of this so urce is quote d in SIMBAD to be z = 1.223 
(fro m Ref. ([S PSS Dat a release T) [2009 ). The Fermi AGN cata- 
log ( |Abdo et al.( |2010b ) ascribes the redhsift 0.678 to the source. 
The source has a powerlaw-type spectrum in the energy range 
0.3-30 GeV with photon index F ^ 2. Significant suppression of 
the source flux above ~ 30 GeV is expected taking into account 
the high redshift of the source. In this respect, the estimate of 
the source flux in the 100-200 GeV band, based on the photon 
detected by the LAT from the source direction, is somewhat sur- 
prising, because it Hes at the extrapolation of the powerlaw from 
the lower energies and is above the expected flux level taking 
into account the expected EBL suppression. At the same time, 
the inconsistency is still very mild due to the large uncertainty 
of the Fermi flux estimate. 



A relatively low level of the source flux makes its detection 
in the ^ > 100 GeV band with MAGIC and CTA difficult if not 
impossible, if the redshift is z = 1.223 (it is, however, detectable 
if the redshift is z = 0.678, as in the Fermi AGN catalog ( |Abdo et| 
[aLl[2Q10 b)). At the same time, a dedicated low-energy-threshold 
ground-based y-ray telescope, like 5@5 would be able to study 
the details of the EBL eff'ects on the source spectrum, because 
the source would be detectable from an energy ~10 GeV. 
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Fig. 10. Lightcurve (top) and spectrum of PKS 195 8- 179 at the 
redshift z = 0.65. Notations are the same as in Fig.[T] 



4. Discussion 

The detection of VHE y-ray emission at energies above 100 GeV 
from blazars with redshifts beyond z = 0.5 carries important im- 
plications for the possibility of studying the cosmological evolu- 
tion of the blazar population and the overall energy output from 
the galaxies, in the form of the EBL. 

Fermi, being able to discover the VHE y-ray signal from 
the high-redshift sources, does not have a capability of study- 
ing the details of the spectral characteristics of these high- 
redshift sources. At the same time, Fermi's detections of these 
sources provides a clear indication for the selection of the 
high-redshift targets for observations with existing and next- 
generation ground-based y-ray telescopes. 

In the previous section we have compared the Fermi flux 
measurements in the £" > 100 GeV band with the sensitivity of 
the ground-based y-ray telescopes. The results from such com- 
parisons indicate that the flux levels of sources listed in table [T] 
are already at or below the sensitivity limit of the current gener- 
ation ground-based telescopes and are, in fact just at the sen- 
sitivity limit of the next generation facility CTA. This means 
that Table [T] provides a more-or-less exhaustive list of the high- 
redshift sources accessible for detection using ground-based y- 
ray telescopes. 

The situation with the limited capabilities for the detection 
of high-redshift sources from the ground could, nevertheless, 
change if a telescope system specially optimized for the reduc- 
tion of the low-energy threshold, like 5@5 would be realized. 
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redshift z= 1.139. Notations are the same as in Fig.[T] 



the redshift z = 1.1. Notations are the same as in Figs. [T] and [3] 



In this case most of the sources listed in Table [T] would be de- 
tectable by 5@5 already at energies of about 10 GeV, with very 
high signal statistics. This would allow a high-quality study of 
the details of the spectral characteristics of the VHE y-ray emis- 
sion from high redshift sources. Lowering the energy threshold 
with 5@5 would open the possibility for detecting a much larger 
number of the high-redshift sources, compared to the just thir- 
teen sources listed in Table [T] 

Low statistics of the VHE y-ray signal from the Fermi high 
redshift blazars does not allow a study of the effect of absorp- 
tion of y-rays through their interactions with EBL photons on a 
source-by-source basis. However, the significance of this eff'ect 
for diff'erent models of cosmological evolution of EBL could 
be evaluated collectively for all thirteen sources from Table JTl 
following the method discussed by |Abdo et al. ( 201Q| ). Fig. [14] 
shows the energies of VHE y-rays from distant blazars as a func- 
tion of the source redshift. The data from the previous analysis at 
lower energies, reported by |Abdo et al.| ( |2Q10 ), are shown in grey 
in the same figure for comparison. The black/grey curves show 
the energies at which the optical depth with respect to pair pro- 
duction on the EBL is T = 1 and r = 3, respectively. This implies 
a flux suppression by a factor of ^ 3 and ^ 20. A naive expec- 
tation is that there should be no VHE photons from the sources 
in the upper part of the plot, much beyond the corresponding 
T = 1 or T = 3 curves. From Fig. 14 one can see that this trend 
of decreasing photon counts beyond the r = 1 curve holds in 
the redshift range z = 0.5 - 1, with only one ph oton beyond the 
solid black curve corresponding to r = 1 for the|Franceschini et 



[aL] ( |2008| ) model of EBL evolution. At the same time, the trend 
seems to be broken for redshifts of about z > 1, where a large 
number of photons beyond the r = 1 curve appears. This large 
number of photons beyond the r = 1 curve is definitely not due 
to the larger source number in this redshift range. Table [T] con- 
tains seven sources at z < 1 and six sources at z > 1. Thus, the 
larger statistics of photons beyond the r = 1 curve should be due 
to a diff'erent reason. 

One possibility is that the determination of redshifts of the 
z > 1 sources listed in the table [T] are not reliable and, in fact, 
the sources are at lower redshifts. In the previous section it is 
mentioned that this might be the case for B2 0912+29. Another 
possibility is that at least one or two photons in the z > 1 part 
of the diagram might still be from the background and not from 
the high-redshift sources. This might be the case for the highest 
energy photon beyond r = 3 curves, which came from the direc- 
tion of RGB J0250+172. This photon has a 0.8% probability of 
being from the background. 

If none of these possible explanations holds (this needs to be 
checked with more data and verification of the redshift measure- 
ments), then the over- abundance of the VHE photons beyond 
the T = 1 curves in the redshift range z > 1 in Fig. [14] should 
be due a physical eff'ect. One possibility is that the EBL level 
at high redshifts is lower, so that the Universe is more trans- 
parent to the VHE photons, than what is implied by the cur- 
rently existing models. Otherwise, new propagation eff'ects, such 
as origin of the y-ray emission in the ultra-high-energy cosmic 
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Fig. 13. Lightcurve (top) and spectrum of PKS 1130+008 at the 
redshift z = 1.223. Notations are the same as in Figs. [T] and |3] 



ray induced cascade (|Essey & Kusen ko[ |2009[ [Murase et al.[ 
'2012! "Aharonian et al., 2012) or even by new physics, like ax- 
ions (Mirizzi et al.||2007) , should be considered. 

In the scenario of |Essey & Kusenko ( 2009| ), UHE protons 
with energies in the EeV range propagate cosmological distances 
and lose energy primarily through the proton pair production 
process. In this case secondary TeV gamma-rays are produced 
by such protons at distances 100-300 Mpc and easily reach the 
Earth. For this model to be valid, one needs a relatively low 
extra-galactic magnetic field with values 1-10 pG everywhere 
on the way of the protons, not only in the voids of large scale 
structure. Apparently, a support for such a scenario is found in a 
recent work of ' Landt| ( |2012 ), in which absorption lines at the 
redshift beyond z = I were found in the spectra of two BL 
Lacs, PKS 0447-439 and PMN J0630-24. One of these BL Lacs 



has been recently observed in the VHE band by HESS (Zech 



et al. 2011). Explanation of the detection of PKS 0447-439 in 



the VHE band within conventional models would be challenging 
if the source redshift is indeed z > I ( Aharonian et al. 1^2012). 
Thorough verification of the redshift measurement for both PKS 
0447-439 and for the z > 1 sources listed in the Table [T] is ex- 
tremely important for the clarificaiton of consistency / inconsis- 
tency of the VHE band detections of the high-redshift sources 
with the current understanding of the mechanisms of formation 
of the VHE y-ray spectra of the sources and of the cosmological 
evolution of the EBL. 



100 




Fig. 14. Energies of VHE photons from BL Lacs and FSRQ as 
a function of redshift. Curves correspond to the optical depth 
T = 1 (black) and r = 3 (grey) with re spect to absorption on 
the EBL. Solid l ine: EBL from Franceschini et al.| ( [2008) ; dotted 
line: EBL from ~ 



Gilmore et al. ( 2012| ); short-dashed line: EBL 
( 2010| ); long-dashed line: EBL from Kneiske et 



from'Finke et al. 
[al. (2004) 
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